1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.02 1.07
## dh0 1.04 1.15
## dh1 1.16 1.55
## dh2 1.14 1.50
## dh3 1.01 1.01
## dl0 1.03 1.13
## dl1 1.05 1.19
##
## Multivariate psrf
##
## 1.16
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1211.32188 | 22616.7719 |
| DIC3 | 1153.83458 | 22718.4531 |
| PWAIC | 42.23985 | 259.0377 |
| WAIC | 1180.55174 | 22741.2220 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.01 1.05
## LDevsum 1.00 1.00
## dh0 1.04 1.15
## dh1 2.70 10.40
## dh2 1.51 2.48
## dh3 2.93 11.81
## dl0 1.09 1.34
## dl1 1.09 1.35
##
## Multivariate psrf
##
## 2.59
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1388.12865 | 22744.0951 |
| DIC3 | 1250.68953 | 22833.8895 |
| PWAIC | 91.90702 | 329.4524 |
| WAIC | 1320.22046 | 22870.8368 |